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ABSTRACT 

Aims. We announce the discovery of a quasar behind the disk of M 31, which was previously classified as a remarkable nova in our 
neighbour galaxy. It is shown here to be a quasar with a single strong flare where the UV flux has increased by a factor of ~ 20. The 
present paper is primarily aimed at the remarkable outburst of J004457+4123 (Sharov21), with the first part focussed on the optical 
spectroscopy and the improvement in the photometric database. 

Methods. We exploited the archives of photographic plates and CCD observations from 15 wide-field telescopes and performed 
targetted new observations. In the second part, we try to fit the flare by models of (1) gravitational microlensing due to a star in M 31 
and (2) a tidal disruption event (TDE) of a star close to the supermassive black hole of the quasar. 

Results. Both the optical spectrum and the broad band spectral energy distribution of Sharov 21 are shown to be very similar to that 
of normal, radio-quiet type 1 quasars. We present photometric data covering more than a century and resulting in a long-term light 
curve that is densely sampled over the past five decades. The variability of the quasar is characterized by a ground state with typical 
fluctuation amplitudes of ~ 0.2 mag around B ~ 20.5, superimposed by a singular flare of ~ 2yr duration (observer frame) with the 
maximum at 1992.81. The total energy in the flare is at least three orders of magnitudes higher than the radiated energy of the most 
luminous supernovae, provided that it comes from an intrinsic process and the energy is radiated isotropically. The profile of the flare 
light curve is asymmetric showing in particular a sudden increase before the maximum, whereas the decreasing part can be roughly 
approximated by a r 5/3 power law. Both properties appear to support the standard TDE scenario where a ~ 10 M Q giant star was 
shredded in the tidal field of a ~ 2 . . . 5 10 8 M Q black hole. The short fallback time derived from the observed light curve requires an 
ultra-close encounter where the pericentre of the stellar orbit is deep within the tidal disruption radius. This simple model neglects, 
however, the influence of the massive accretion disk, as well as general-relativistic effects on the orbit of the tidal debris. Gravitational 
microlensing probably provides an alternative explanation, although the probability of such a high amplification event is very low. 

Key words. Quasars: general - Quasars: individual: 1004457+4123 - Galaxies: individual: M 3 1 - Gravitational lensing - Black hole 
physics 
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03 | 1. Introduction of the source, as is the case for quasars, takes a long time. 

This also holds for recurrent events, e.g. novae, with intrinsi- 

Temporal variability is one of the most conspicuous properties call y shorter time scales where h is necessary to cover also the 

for several classes of interesting astrophysical objects. Owing to wlde 8 a P s between the single events. Presently, data mining in 

the unprecedented combination of sky coverage and photometric archives, most notably in the plate archives from large Schmidt 

accuracy, discoveries from the Large Synoptic Survey Telescope telescopes, remains the only approach if the light curves have 

(LSST), the Panoramic Survey Telescope & Rapid Response t0 cover a tlme lnterval of decades in the rest frame, at least for 

System (Pan-STARRS), the Palomar-QUEST (PQ) survey, or high-redshift quasars. 

the Palomar Transient Factory (PTF) will provide great advances In the context of the search for and the identification of 

in the understanding of variable processes, especially of rare optically variable star-like sources in the field of the bright 

transient phenomena (Gezari et al. I2008I Strubbe & Quataert Local Group spiral galaxy M31 (Pietsch et al. 2005a; Henze 

2009; Quimby et al. 2009). However, given the nature of the et al. 2008 ), our interest was pointed toward the apparent nova 

problem, the creation of the observational database for investi- J004457+4123, originally discovered by Nedialkov et al. (1996 ) 

gating the variability on long time scales in the reference frame and described in more detail by Sharov et al. ([1998 ). The light 
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curve presented by these authors clearly shows a strong bump 
with a maximum brightening by more than 3 mag in the year 
1992 while the source remained constant both in the 23 years 
before and the 5 years after. Sharov et al. suggest that it is a 
"remarkable nova" in M31, but underscore that it "differs dra- 
matically from typical representatives of this class of objects". 
Following the terminology of these authors (nova 21), we de- 
note the object as Sharov 21 throughout this paper. A possible 
X-ray counterpart was first discussed by Pietsch et al. (12005b). 
who searched for supersoft X-ray counterparts of optical novae 
in M 31 and identified Sharov21 with the source [PFH2005] 601 
of their catalogue of XMM-Newton EPIC X-ray sources (Pietsch 
et al. I2005al ) and with the hard R OSAT source [SHL2001] 306 
from the catalogue of Supper et al. (2001 ). Based on the hardness 
of the X-ray source and the peculiar optical light curve, Pietsch 
et al. d2005bl > speculate that Sharov 21 "may not be a nova at 
all". Here we present, for the first time, optical follow-up spec- 
troscopy which reveals Sharov 21 to be a quasar. 

From the very beginning of the investigation of active galac- 
tic nuclei (AGN), variability is known to be a diagnostic property 
of this object class and has been successfully used as a crite- 
rion for the selection of quasar candidates in a number of stud- 
ies (e.g., Kron & Chiu 119811 Majewski et al. I199U Hawkins 
& Ve ron [19931 Meusinger et al. [20021 [20031 Rengstorf et al. 
2004). AGN originally misclassified as variable stars are nei- 
ther unprecedented nor unexpected. The most famous case is the 
prototypical blazar BL Lac, discovered by Cuno Hoffmeister in 
1930. However, the misclassification of a luminous quasar as a 
nova is highly remarkable because it indicates a singular, strong 
outburst which points toward a rare and interesting transient phe- 
nomenon. 

It has long been understood that the observed flux varia- 
tions of AGNs hold keys to the structure of the radiation source. 
The physical mechanisms behind these fluctuations are however 
still poorly understood. Frequently discussed scenarios for the 
origin of the observed optical/UV broad-band long-term (non- 
blazar) variability related to massive or supermassive black holes 
in galaxy centres include various processes such as instabilities 
and non-linear oscillations of the accretion disk (Taam & Lin 
119841 Abramovicz et al. 119891 Honma et al. I199U Kawaguchi 
et al. 1998), multiple supernovae in the starburst environment 
(Terlevich et al. 119921 Cid Fernandes et al. 119971 ). microlens- 
ing of the accretion disk or the broad line region by com- 
pact foreground objects (Chang & Refsdal, 119791 Irwin et al. 
[19891 Hawkins [T9911 Schneider IT9931 Lewis & Irwin [19961 
Zackrisson et al. 2003 ), the disruption of a star which passes 
within the tidal radius of the supermassive black hole (Hills 
[19751 Rees [19881 [19901 Komo ssa & Bade [19991 Komossa & 
Meritt 2008; Gezari et al. 120081 ). star-star collisions in the dense 
circumnuclear environment (Torricelli-Ciamponi et al. [2000), 
and interactions of the components in a supermassive binary 
black hole (Sillan paa et al. [T9881 Lehto & Valtonen [19961 Katz 
[19971 Liu & Chen[2007]). 

The present paper is aimed at the highly peculiar light curve 
of the quasar Sharov 21 which is worth detailed investigation. 
We present the optical spectrum and a significantly improved 
light curve and discuss possible scenarios for the strong out- 
burst. The observations are described in Sect. [2] The spectrum 
and other basic properties are analysed in Sect.0 The outburst is 
the subject of Sect.[4] Two models are discussed in detail: gravi- 
tational microlensing and a stellar tidal disruption event; alterna- 
tive scenarios are briefly summarized as well. Section[5]gives the 
conclusions. Standard cosmological parameters Hq = 71 km s _I 
Mpc~', Q. m = 0.27, Qa = 0.73 are used throughout the paper. 



2. Observational data 

2.1. Spectroscopy 

The optical spectrum was obtained with the Double Imaging 
Spectrograph (DIS) on the 3.5-m telescope at Apache Point 
Observatory (APO) in New Mexico, USA during a campaign 
to follow-up X-ray sources in M 3 1 . Two exposures were taken 
for a total of 4 500 seconds. For the blue spectral range (3 200 to 

5 500 A), the B 400 reflectance grating was used with a disper- 
sion of 1.85 A per pixel yielding a nominal resolution of about 
7 A in combination with a 1'.'5 entrance slit. The R 300 grating, 
with a dispersion of 2.26 A per pixel, gives a resolution of 8 A 
for the red part (5 000 to 10 100 A). The spectra were taken at 
UT 0300 on 2007-1 1-09. The observing conditions were excel- 
lent through the night. 

The spectra were reduced, wavelength calibrated, and flux 
calibrated using the standard IRAF routines (ccdproc, iden- 
tify, sunsfunc, apextract, apall). Wavelength calibration was per- 
formed using HeNeAr lamp exposures taken just before the ob- 
ject exposures, and flux calibration was performed using a spec- 
trum of the spectrophotometric standard star BD+28-421 1. 

2.2. Optical photometry: long-term light curve 

The light curve published by Sharov et al. ( 119981 ) is based on 
B band observations taken with four telescopes between 1969.0 
and 1997.7 with a good coverage of the outburst phase. The 
present study is aimed at an extended and better sampled long- 
term light curve. We exploited several data archives and com- 
bined the results with the data available in the literature. In ad- 
dition, targetted new observations for another 16 epochs in the 
years 2006 to 2009 were taken with the CCD Schmidt camera of 
the Tautenburg 2 m telescope and with the focal reducer camera 
CAFOS at the 2.2 m telescope on Calar AltcQ Spain. Most of the 
archival photographic plates were digitized in the frame of the 
present work using the Tautenburg Plate Scanner (Brunzendorf 

6 Meusinger, 1999 ) for the Tautenburg Schmidt plates, the high- 
quality commercial scanner at the Asiago observatory (Barbieri 
et al. 2003 ) for the Asiago plates, and the Microtek ScanMaker 
9800XL for the Sonneberg astrograph plates. The Calar Alto 
plates were scanned for the Heidelberg Digitized Astronomical 
Plates (HDAP) project using a Heidelberger Druckmaschinen 
Nexscan F4100 professional scanner and are available from the 
German Astrophysical Virtual Observatory (GAVOfl 

A summary of all used observations from the last six decades 
is given in TableQ] (CFHT = Canada France Hawaii Telescope, 
INT = Isaac Newton Telescope, WFS = wide-field survey). N t 
is the total number of all single exposures, N e the number of 
epochs in the light curve where the quasar has been measured. 
The last column gives the source of the photometric reduction: 
(1) this work, (2) Vilard ell et al. d2006b . (3 ) Mon et al al. d2003l . 
(4) Massey et al. J20061 ). (5) Sharov et al. (TT9981 ). 

Altogether, the light curve data pool contains more than 
1 100 single observations from 15 telescopes. Included are the 
B magnitudes published by Sharov et al. ( 1998 ) for 84 epochs, 
by the Local Group Galaxies Survey (LGGS; Massey et al. 
2006) for 2 epochs, and by Vilardell et al. (2006) binned here 
into 6 epochs. For the other observations, the photometric re- 



1 The Calar Alto Observatory of the Centro Astronomico Hispano 
Aleman, Almeria, Spain, is operated jointly by the Max-Planck-Institut 
fur Astronomie and the Instituto de Astrofisica de Andalucfa (CSIC). 

2 http ://dc . zah .uni-heidelberg . de 
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Table 1. Observational material for the construction of the light 
curve. 



telescope 




N e 


years 




source 


(a) Digitized photographic plates 










Asiago Schmidt 


24 


8 


1968.8- 


-1993.1 


(1) 


Calar Alto Schmidt 


43 


15 


1983.0- 


-2000.7 


(1) 


Calar Alto 1.2 m 


8 


5 


1976.7- 


-1982.6 


(1) 


Palomfir Schmidt 


4 


4 


1948.7- 


-1989.7 


Q).(3) 


Sonneberg 40 cm 


9 


1 


1992.2- 


-1992.6 


(1) 


Tautenburg Schmidt 


362 


77 


1961.5- 


-1997.0 


(1) 


(b) CCD observations: 












Calar Alto 2.2 m 


3 


2 


2008.7- 


-2009.3 


(1) 


CFHT3.6m 


1 


1 


1993.8 




(1) 


INT (WFS) 


5 


1 


1998.8 




(1) 


INT 


522 


6 


1999.7- 


-2003.7 


(2) 


KittPeak4m 


10 


2 


2000.8,2001.7 


(4) 


Skinakas 60 cm 


5 


1 


2007.6 




(1) 


Tautenburg Schmidt 


30 


14 


2006.1- 


-2009.7 


(1) 


(c) Original data from Sharov et al. (1998 ): 






four other telescopes ~ 


150 


84 


1969.0- 


-1997.7 


(5) 



duction was done in the frame of the present study. We used 
the Source Extractor package (Bertin 1 19961 1 for object selec- 
tion, background correction, and relative photometry, and the 
LGGS catalogue for the photometric calibration. The reduction 
was performed under ESO MIDAS. Because of the strongly in- 
homogeneous background across the disk of M 3 1 (Henze et al. 
2008), the photometric calibration was done locally on a 8' X 8' 
subimage around Sharov 21 where typically ~ 100 + 50 calibra- 
tion stars from the LGGS were identified. Note that the magni- 
tudes given by Sharov et al. ( 1998) have also been derived from 
standard stars close to the target. The blue magnitude for the 
Palomar POSS 1 plate (1953-1 0-09) is taken from the USNO- 
B1.0 catalogue (Monet et al. 120031 1 whereas the magnitudes 
given there for the POSS 2 plate do not agree with the impres- 
sion from the visual inspection of the images. We re-reduced 
the image cutouts from the Digitized Sky Survey (DSS2) and 
derived B = 20.73 + 0.23 and R = 19.66 + 0.25. An early 
deep plate taken with the 1.2 m Samuel Oschin Telescope on 
1948-09-29 is shown in the Hubble Atlas of Galaxies (Sandage 
119611 1. The visual inspection shows that Sharov 21 is detected 
with B = 20.3 + 0.3. 

With very few exceptions the observations were made 
through filters reproducing the Johnson UBVR system. About 
80% of the data points in the light curve are from observa- 
tions in the B band. From several pairs of observations, taken 
at nearly the same epoch but with different filters, we derive re- 
lations between colour indices and the B band magnitude. The 
results (Fig. [3) are used to obtain B from observations in the 
other bands. Because quasar variability is usually not achro- 
matic (see Sect.O, such colour relations are more suitable than 
single-epoch colour indices. Photographic magnitudes m pg from 
blue sensitive emulsions without filter were transformed by B ~ 

TOpg +0.1. 

There are two basic limitations with regard to the final set of 
data. First, with B > 20 for most of the time, Sharov 21 is too 
faint to be detected in every plate archive. A high fraction of ob- 
servations yields therefore upper limits only. In other cases the 
detections are close to the plate limit resulting in relatively large 
photometric errors. Second, the epochs of the observations are 



not regularly distributed but show strong clustering. We stacked 
the images taken with the same telescope within typically sev- 
eral days to a few weeks applying a quality-weighting factor 
(Froebrich & Meusinger 2000) to obtain deeper images with 
improved signal-to-noise ratio. This procedure was not applied 
however for the outburst phase in 1992 where we are interested 
in a high temporal resolution. 

The final light curve (Fig.Q} comprises magnitudes at 221 
detection epochs but still suffers from several gaps. In particu- 
lar, no data are available for the early rising phase of the out- 
burst between March and August 1992. We checked the Wide 
Field Plate Datebas43but found no entries for this time. Also the 
search in the plate archives of the Baldone Schmidt telescope 
(Alksnis et al. [19981 and of the 100/300 cm Schmidt telescope 
of the Kvistaberg Observatory revealed no observations of M 3 1 
during that time. 



Table 2. Detection limits on plates taken before 1950. 



telescope 


plate 


date 


"^pg,lim 






yyyy-mm-dd 




Bruce astrograph LHK 


23a 


1900-09-14 


18.8 


Bruce astrograph LHK 


265a 


1901-08-18 


18.0 


Yerkes 24 inch reflector 


? 


1901-09-18 


19.5 


Bruce astrograph LHK 


649a 


1903-01-15 


18.0 


Bruce astrograph LHK 


842a 


1903-09-27 


18.0 


Bruce astrograph LHK 


1384a 


1905-12-26 


18.0 


Waltz reflector LHK 


198 


1907-11-02 


17.5 


Waltz reflector LHK 


369 


1908-08-20 


18.3 


Waltz reflector LHK 


603 


1909-10-19 


18.5 


Waltz reflector LHK 


4484 


1934-09-17 


18.5 


Bruce astrograph LHK 


7163a 


1949-09-20 


18.2 



It is useful to check also older historical observations of 
the Sharov 21 field which are not deep enough to detect the 
quasar in its faint stage but would allow discovering a previ- 
ous outburst. Table|2]lists those observations from the 40/200 cm 
Bruce double-astrograph and the 72 cm Waltz reflector of the 
Landessternwarte Heidelberg-Konigstuhl (LHK) with m pg; ij m ~ 
18. A reproduction of a plate taken in 1901 with the 24 inch re- 
flector at Yerkes observatory is shown by Hubble (1929); from 
the visual inspection we estimate a detection threshold m ~ 19.5. 
Each of these observations excludes the occurrence of a flare 
similar to that of 1992 for at least several tens of days around 
their dates of exposure. 



3. General properties of Sharov 21 

In Sect. l3.4l below we demonstrate that the previous classifica- 
tion of Sharov 21 as a remarkable nova in M31 has to be re- 
jected. Our optical spectrum, presented below, clearly reveals 
the source to be a quasar. The most important properties of this 
quasar are summarized in Table[3] fi/2 is the time interval for the 
decline from maximum flux to half the maximum. Remarks: (1) 
position from LGGS, (2) ground state/maximum, (3) observer 
frame/quasar rest frame, (4) based on C iv line, (5) mean value 
for the ground state. 



3 Both the WFPDB version available via VizieR, (Tsvetkov et al. 
1 19971 1, and the updated version available via the search browser devel- 
oped in Sofia (http://draco.skyarchive.org/search) were used. 
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Fig. 1. Long-term B light curve from the data summarized in TableQ] (symbols plus error bars; no photometric errors available for 
the data from Sharov et al. Q998 )). Dotted vertical lines: upper limits from images on which the object is not detected. For lucidity, 
only a small fraction of the upper limit data is shown. 



Table 3. Basic properties of Sharov 21 (remarks: see text). 



Measured and derived quantities 




Remark 


RA (2000) 


00 h 44 m 57 s 94 


(1) 


Dec (2000) 


+41°23'43'.'9 


(1) 


redshift z 


2.109 




projected distance from M 31 centre 


26' 




apparent magnitude B 


20.5/17.2 


(2) 


foreground dust reddening E(B — V) 


0.2 mag 




absolute magnitude M B 


-27.5/ - 30.7 


(2) 


date of the maximum (year/JD) 


1992.81/2448918 




ti/2 decline (days) 


15/5 


(3) 


black hole mass M^h 


5 10 s M G 


(4) 


log (AW erg s" 1 ) 


46.6 


(5) 


Eddington ratio Aoi/Add 


0.60 


(5) 



3. 1 . Optical light curve and general remarks on variability 

The light curve from the data discussed in Sect. [2] is shown in 
Fig.Q] Sharov et al. ( 1998) note that the quasar was "nearly con- 
stant from 1969 through 1991 with B « 20.5 and returned to 
this value one or two years after the outburst". Compared with 
the original data from Sharov et al., we (1) basically confirm 
their finding, (2) extend the covered time interval by about one 
decade in each time direction, and (3) fill some broad gaps (e.g., 
between the years 1984 and 1990 and between the beginning of 
1993 and 1994). Based on the better sampling of our data, in- 
cluding the upper limits, the possibility of outbursts in intervals 
that were not covered by observations is thus significantly re- 
duced. Hence, we conclude that the light curve can be divided 
into (a) the faint state (with B = 20.52), which can be con- 
sidered as the ground state, and (b) a single outburst, or flare, 
lasting ~ 2 yr (JD ~ 2 448 500 ... 2 449 300) where the quasar 
was 3.3 mag brighter in the maximum. The flare (Fig.[2j shows 
a slightly asymmetric profile with three phases: (1) a gradual in- 
crease between JD -2448 500 and 2448 880 with a gap in the 
light curve between March and August 1992, followed by (2) an 
abrupt rise to the maximum at JD ~2 448 918, and (3) a quasi- 
exponential decline to the ground state at JD ~2 449 300. The 
interpretation of the outburst will be the subject of Sect. [4] 




4BBOO 49000 
J. D. - 2400000 



Fig. 2. The light curve in the outburst. Open symbols: uncertain 
data, dashed lines: upper limits, dotted horizontal line: ground 
state. 



To evaluate the variability in the ground state we com- 
pare Sharov 21 with quasars from the Tautenburg-Calar Alto 
Variability and Proper Motion Survey (VPMS; Meusinger et al. 
2002, 2003). For Sharov 21 we have a B standard deviation <x B = 
0.27 mag from the Tautenburg data (0.26 mag from all data). For 
the VPMS quasars with similar redshifts (z = 2.1 +0.2) and com- 
parable (extinction-corrected) mean magnitudes (B = 19.7 + 0.2) 
we have cr B = 0.26 mag in the VPMS field around M3 (8 
quasars) and 0.29 mag in the field around M92 (6 quasars). We 
conclude that the flux variability of Sharov 21, in its ground 
state, is not unusually strong. 

Variations of the B band flux of Sharov 2 1 are correlated with 
colour changes. The observed relations (Fig. 0) are qualitatively 
in agreement with the typical properties of quasars. A harden- 
ing of the optical/UV continuum during the bright phase is in- 
dicated by multi-frequency monitoring of selected AGNs (Cutri 
et al. [19831 Edelson et al. [19901 Paltani & Courvoisier [19941 ) as 
well as by statistical studies of AGN ensemble variability (Di 
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17.5 18 18. 5 19 19.5 20 20.5 21 

B 

Fig. 3. Colour indices of Sharov21 as a function of its B mag- 
nitude from quasi-simultaneous observations in both bands. 

Clemente et al. [19961 Cri stiani et al. [T997l Trevese et al. l200Tl 
Vanden Berk et al. 2004). This trend has been confirmed also by 
multi-epoch spectroscopy of quasars from the Sloan Digital Sky 
Survey (SDSS) where it was shown that the emission lines are 
considerably less variable than the continuum, being stronger in 
the faint stage, relative to the continuum, than in the bright phase 
(Wilhite et al. 2005 ). For Sharov21, the U band is dominated by 
the strong Lyman or/N v line (Fig.[5]l. The contribution of the C iv 
line to the flux in the B band is much smaller, and the other bands 
are nearly pure continuum. The colour indices B - V and B - R 
are hence expected to become bluer when the quasar becomes 
brighter, while U - B becomes redder at the same time. Wilhite 
et al. present the colour differences between the bright and the 
faint phase as a function of redshift (their Fig. 14) indicating that 
A(u - g) has a local minimum at z ~ 2 while A(g - r) has a peak 
at the same redshift. Sesar et al. (2007) found a similar result 
from the multi-epoch photometric data of quasars in the SDSS 
stripe 82 with respect to the g and r bands. In an ongoing study 
of quasar variability in the SDSS stripe 82 (Meusinger et al., in 
preparation) we confirmed the relation shown by Sesar et al. and 
derived corresponding relations for the other bands which are in 
line with the results from Wilhite et al. (2005 ). 

3.2. Position behind M 31 

Figure|4]reveals that Sharov21 is seen through the disk of M 31. 
The ellipses represent the D25 isophotes of M31, M32, and 
NGC205, respectively, according to the RC3 (de Vaucouleurs 
et al. fl99lT ). For a distance of 750 kpc for M 31 (Vilardell et al 



2006 1 their table 1), an inclination angle i = 77°, and a posi- 
tion angle of the major axis of 35° from the RC3, its projected 
distance from the centre of 26' (5.7 kpc) corresponds to a galac- 
tocentric distance of R = 16 kpc in the midplane of M31. A 
radius of ~ 30 kpc is a realistic assumption for the extent of the 
bright disk (e.g. Racine 1991; Ferguson et al. 2002; Irwin et al. 
2005 ). The "Catalogue of Quasars and Active Galactic Nuclei 
(12th Ed.)" (Veron-Cetty & Veron 2006) lists only three quasars 
within 100' from the centre of M31. However, Sharov21 is, to 
our knowledge, the first quasar detected behind the disk of M 3 1 . 




12 11 10 9 



RA ( 200:0 ) [deg] 

Fig. 4. Map of the known quasars in the M31 field. Open 
squares: quasars from Veron-Cetty & Veron (120061 ). framed as- 
terisk: Sharov21. 



3.3. Foreground reddening 

A foreground spiral galaxy is expected to produce a substan- 
tial reddening of a background quasar (Ostman et al. 2006). The 
"Galactic Dust Extinction Service'Qof the NASA/IPAC Infrared 
Science Archive, which is based on the method pioneered by 
Schlegel et al. ( fl998l ), provides E(B - V) = 0.20 mag for the 
position of Sharov21. Individual reddening values for a large 
set of globular clusters in M31 were derived by Barmby et al. 
(2000) and Fan et al. (120081 yielding E(B-V) = 0.12 + 0.03 mag 
for the six clusters within 6' from Sharov21. The reddening of 
the quasar must be stronger since only a fraction of the clus- 
ters is expected to be located behind the disk of M 3 1 . Finally, 
a simple model for the radial dependence of the extinction in 
M31 derived by Hatano et al d 19971 ) yields A B = 0.68; i.e., 
E(B - V) = 0.17 mag for the standard Milky Way extinction 
curve (Savage & Mathis [T979l l, which seems to be valid also for 
M31 (Barmby [2000]!. Here we adopt E(B - V) = 0.20 mag for 
the total foreground reddening of Sharov21. 

3.4. Optical spectrum 

The optical spectrum is shown in Fig. [5] For comparison the 
composite spectrum of "normal" quasars from the Sloan Digital 

4 http://irsa.ipac.caltech.edu/applocations/DUST/ 
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Sky Survey (SDSS; Vanden Berk et al. l200Tb is plotted (thin 
smooth curve) shifted to the redshift of Sharov21, the spectra 
are normalized at A 4500 A. We derived a redshift of z = 2.109 
both from the fit of the SDSS composite and directly from the 
wavelengths of the narrow components of the Lyman a and N v 
lines. Compared with the SDSS composite Sharov21 has a red- 
der continuum. We de-reddened the spectrum for foreground 
(z = 0) extinction adopting the Milky Way extinction curve. 
Good agreement with the mean SDSS quasar spectrum is found 
for E(B-V) = 0.2 mag (Fig. [5] bottom), which is perfectly in line 
with the reddening value from the NASA/IPAC Infrared Science 
Archive (see above). 

The de -reddened spectrum of Sharov21 is that of a typi- 
cal type 1 quasar. There is no evidence of unusual spectral fea- 
tures indicating a peculiar nature of Sharov21. Compared with 
the SDSS composite the reddening corrected spectrum shows 
a stronger Fe bump at A ~ 7300. . .8300 A (observer frame) 
which points towards a relatively high Eddington ratio (Dong et 
al. 120091 . The Cm] A 1909 A line appears slightly weaker, but 
note that the line coincides with the Na i foreground absorption 
at AA 5890, 5896 A. We notice further a weak unidentified ab- 
sorption line at the position of the N iV/i 1486 A line. 




4000 5000 6000 7000 B000 



uaue length A 

Fig. 5. Observed (top) and foreground extinction-corrected (bot- 
tom) optical spectrum of Sharov21 (observer frame), not cor- 
rected for telluric absorption. 



3.5. Other wavelength regimes 

M31 has been observed at radio wavelengths both as part of 
larger surveys and as the focus of dedicated programmes (see 
e.g., Gelfand et al. 2004; their Table 1). In no case, a radio coun- 
terpart was detected at the position of Sharov21. With a flux 
density limit of 0.15 mJy, the VLA survey by Braun ( 1 19901) is 
the deepest one at 20 cm continuum ( 1 .465 GHz). At lower radio 
frequencies a deep survey of the M 3 1 field was performed by 
Gelfand et al. d2"U04l 120031 ) with the VLA yielding a flux den- 
sity limit of 3 mJy at 325 MHz. Within d ~ 2' from Sharov21, 
the 325 MHz catalogue lists the sources GLG043, 045, 050, 
and 051. However, with d > 1' all four radio sources are ob- 



viously not related to Sharov21. The non-detection in the Braun 
survey implies an upper limit for the radio-loudness parameter 
R, i.e. the ratio of the 5 GHz radio flux density to the 2500 
A optical flux density in the quasar rest frame (Stocke et al. 
119921 ), log/?* < 0.7 for a radio spectral index ctr = -0.3 and 
log/?* < 0.9 for a R = -0.5. With t he thre shold log/?* > 1 
for radio-loud AGNs (e.g., White et al. 2000} Sharov21 is not a 
radio-loud flat-spectrum quasar. 

We re-analysed the archival XMM-Newton and ROSAT 
data. For computing fluxes from instrument dependent count 
rates, we used an absorbed power-law model with a generic pho- 
ton index of 1.7 (see also Pietsch et al. 2005a). We adopted the 
foreground extinction of E(B - V) — 0.20 mag, derived form 
the optical data, which translates to a Nn of 1.1 10 21 crrT 2 fol- 
lowing Predehl et al. ([1995). Based on this spectral model, we 
used the source count rates given by Pietsch et al. (2005aj( for 
[PFH2005] 601 in the XMM-Newton observation 0151580401 
(2003-02-06) to estimate an unabsorbed flux of (5.8 + 1.9) 10" 14 
erg s _1 cm" 2 in the (0.2 - 10.0) keV band and the monochromatic 
flux F v = (4.5 + 1.5) 10~ 9 Jy at v = 1.7 10 17 Hz. 

The ROSAT data of the Sharov 21 fielc0 consists of 12 obser- 
vations with the Position Sensitive Proportional Counter (PSPC) 
and 3 observations with the High Resolution Imager (HRI). The 
data analysis was done under ESO MIDAS within the EXSAS 
context. We performed source detection around the position of 
Sharov 21 on the original event files and computed count rates 
and 3<r upper limits for all observations. There is just one 3<x 
detection of an X-ray source in the data set, identical with 
[SHL2001] 306, which is supplemented by upper limits for the 
rest of the observations. The flux estimated from the detection 
is consistent, within the errors, with the XMM-Newton data. 
Although the ROSAT observations were performed around JD 
2449000, i.e. during the decline of the UV flare, no significant 
X-ray variability is detected. Note however, that the only ROSAT 
detection is a very faint off-axis PSPC detection and due to the 
large positional error circle of this instrument we can not assume 
a doubtless correlation with [PFH2005] 601. 
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Fig. 6. SED of Sharov 21 in the restframe (open symbols with 
downward arrows: upper limits). For comparison the mean SED 
(Elvis et al. 119941 ) is shown for radio-quiet (solid) and radio-loud 
quasars (dotted), normalized at A 1415 A. 



5 see http://www.xray.mpe.mpg.de/cgi-bin/rosat/seq-browser 
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We also checked images from the Deep Imaging Survey with 
the Galaxy Evolution Explorer (GALEX) in the near ultravio- 
let (A eB ~ 2270 A) and in the far ultraviolet (A eS ~ 1520 A), 
but no counterpart could be identified within a radius of ~ 5 ". 
For detections in the near infrared we searched in the 2MASS 
catalogue (Skrutskie et al. 2006 ), again without a clear-cut iden- 
tification of a counterpart. From simple statistics of the sources 
around the position of Sharov 21a flux limit of 1 yuJy is estimated 
for the ultraviolet bands and of 0.25 mJy in the K band. In Fig. [6] 
we compare the extinction-corrected broad band spectral energy 
distribution (SED) of Sharov 21 with the mean SED for normal, 
nonblazar quasars from Elvis et al. d 19941 1. 

3.6. Black hole mass and Eddington ratio 

The black hole mass is estimated from the C iv line width as a 
measure of proxy for the velocity dispersion of the broad emis- 
sion line gas in combination with a radius-luminosity (R-L) re- 
lationship for the emission region. Significant progress has been 
achieved over the last years with the calibration of the R-L rela- 
tion ( Vestergaard [20 021 120091 Corbett et al. [20031 Warne r et al. 
120031 Peterson et al. l2004l Vestergaard & Peterson 12006b . 

We use equation (8) from Vestergaard & Peterson 
(2006) which is based on the line dispersion cr line (Civ) and 
the monochromatic continuum luminosity ALx at 1 350 A. 
Following the method outlined by Peterson et al. (120041) and 
Vestergaard & Peterson (2006) we obtain cr line (Civ) = (2.9 + 
0.5) 10 3 km s _1 . The luminosity of the continuum at A = 1 350 A 
(restframe) is derived from the extinction-corrected mean B 
band flux in the ground state after correction for the contribu- 
tion from the emission lines (~ 12%) and assuming a standard 
power-law continuum F,\ oc A~ (2+a) with a = -0.5. We find 
AL A (1 350 A) = 1.07 10 46 erg s _1 which finally yields a black 
hole mass of M b h ~ 5 10 8 M o . Vestergaard & Peterson (2006) 
give a standard deviation of +0.33 dex for the scaling relation. 
Allowing further for the uncertainties of the line width and of the 
continuum luminosity, the total uncertainty of Mbh is roughly a 
factor of 3. 

Adopting the bolometric correction &b c (l 350 A) = 4 from 
Richards et al. (2006; their Fig. 12), the monochromatic lumi- 
nosity from above corresponds to Lboi = ALx ■ K c = 4.3 10 46 
ergs -1 ~ 10 13 L o and to an Eddington ratio e = Lb i/£ e dd = 
3.3 10 8 (M hh /M Q Y l ~ 0.6 for the ground state. The overwhelm- 
ing majority of SDSS quasars with z ~ 2 radiate at between 10% 
and 100% of their Eddington luminosity (see Vestergaard 2009) 
with a narrow e interval at high luminosities (Kollmeier et al. 
120061 Shen et al. 120081 Onken & Kollmeier [2008). According 
to Onken & Kollmeier, quasars with log (Lb i/erg s ) > 46.5 
and z = 1.1 . . .2.2 are characterized by (log e) = -0.65 + 0.35 
whereas e is smaller and shows a broader distribution for low- 
luminosity AGNs. Gavignaud et al. ( 120081 ) find log e - -0.97 + 
0.28 [log(Lboi/erg s _1 ) - 45] over the interval log(Lboi/i«) = 
45 . . .48.5. As can be seen from their Fig. 4, a value of e ~ 0.6 
is relatively high but not unusual for quasars of the redshift and 
luminosity of Sharov 21. However, in the maximum of the out- 
burst, where the B band flux is a factor of ~ 20 higher, the lu- 
minosity of Sharov 21 corresponds to a highly super-Eddington 
regime. 

It has been questioned for a couple of reasons whether the 
mass estimates from scaling relations over-predict the black hole 
masses by about an order of magnitude. For arguments in favour 
of the mass estimates we refer to the discussion by Vestergaard 
(2004). In addition we note that a significantly lower black 



hole mass would require an unusually high Eddington ratio for 
Sharov 21. 



4. The outburst 

4.1. Sharov 21 -like UV flares are rare 

UV Flux variations by one or a few magnitudes are observed 
both in low-luminosity AGNs like NGC 5548 (Ulrich et al, IT997l 
and references therein) or in blazars (Sect. 14.4) . Contrary to the 
burst of Sharov 21, the B band flux of high redshift radio-quiet 
quasars typically varies by a few tenths of a magnitude as is ob- 
served in the ground state of Sharov 21. The second remarkable 
difference is the fact that the strong activity of our quasar is lim- 
ited to a short time interval of ~ 1 year (observer frame). During 
the remaining ~ 47 yr covered by the light curve (i.e., the ground 
state), the mean flux variation is a factor of ~ 16 smaller than the 
maximum fluctuation in the outburst. In contrast, both strongly 
variable low-luminosity AGNs like NGC 5548 and optically vi- 
olently variable (OVV) quasars show a more or less steady up 
and down variation. 



0.2 0.4 

aCmag) 




log Cr/days) 

Fig. 7. Comparison of the variability of Sharov 21 with the 
VPMS quasars. Top: maximum fluctuation AB max vs. standard 
deviation <x (framed asterisk: Sharov 21; <x refers to the ground 
state). Bottom: single object structure functions for Sharov21 
(solid line at the top) and the VPMS quasars (thin solid lines), 
and VPMS sample-averaged structure function (dotted line). 
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To get an idea how usual or unusual the outburst of Sharov 2 1 
is, we first consider the light curves of the quasars from the 
Variability and Proper Motion Survey (VPMS; Meusinger et al. 
2002, 2003 ). The time baseline of the observations is nearly the 
same (from ~ 1960 to 2008) as for Sharov21, the numbers of ob- 
servations per quasar are however much lower, typically ~ 50. 
The VPMS quasar sample is highly complete down to B ~ 20. 
Among the 321 AGNs in the VPMS there are 10 AGNs with 
maximum fluctuations AB max = B max - B m [ n > 2 mag. For 
the majority of them (80%), the variability is characterized by 
a more or less monotonic variation over the time baseline. Only 
for two AGNs the maximum fluctuations can be attributed to 
burst-like features. One is the strongly variable object CC Boo, 
a Seyfert galaxy at z = 0.17 (Margon & Deutsch [19971 ), the 
other one is a radio-quiet quasar at z = 1.08. However, with 
AB max = 2.46 mag and 2.50 mag, respectively, the burst ampli- 
tudes are considerably smaller than for Sharov 21 and there are 
strong fluctuations also in other parts of the light curves. 

In Fig.|7]we compare variability properties of Sharov 21 with 
those of the VPMS quasars. The top panel shows the maximum 
fluctuation amplitude versus standard deviation of the B magni- 
tudes in the ground state for Sharov 2 1 , compared with 98 VPMS 
quasars of comparable redshifts (z = 2.1 ± 0.5). As we cannot 
reasonably distinguish between a ground state and a higher state 
for the majority of the quasars we simply use the standard devi- 
ation of all data in the light curve as a proxy. (A natural conse- 
quence is the increase of cr with the maximum amplitude.) The 
light curves of the two strongly variably VPMS quasars marked 
by open squares in Fig.|7]are clearly dominated by smooth long- 
term variations over decades. A popular statistical tool for the in- 
vestigation of quasar variability is the first order structure func- 
tion S 2 (t) = {[ m(t + t ) - m(t)] 2 ) t (e.g., Simonetti et al. 119851 
Kawaguchi et al. 1998) where r is the time-lag between two ob- 
servations in the quasar restframe and the angular brackets de- 
note the time-average. The structure function represents a sort 
of running variance (as a function of the time-lag) and contains 
therewith information about the time scales of the involved vari- 
ability processes. The most important conclusion from Fig.|7]is 
that the flare of Sharov 21 is singular and without comparison in 
the long-term variability data of the VPMS quasar sample. 

Excellent data for the statistical study of quasar variabil- 
ity has been provided from stripe 82 of the Sloan Digital Sky 
Survey (SDSS) for ~ 10 4 quasars in five colour bands over 
~ 7 years (e.g., Sesar et al. 2007). Using the SDSS quasar cat- 
alogue (Schneider et al. 120071) we identified 8311 quasars in 
the Light and Motion Curve Catalogue (LMCC; Bramich et al. 
2008 ("from ~ 249 square degrees "of the SDSS stripe 82. No 
high-redshift quasars with z > 2 were found with amplitudes 
in the u and g bands > 1.5 mag. Allowing for the whole red- 
shift range, the two SDSS quasars with the highest amplitudes 
in the g band are SDSS J001 130.0+005751.8 (z = 1.49) and 
SDSS J21 1817.37+001316.8 (z = 0.46) with Ag m . dX = 3.2 and 
2.7 mag, respectively. Both are bright polarized flat-spectrum ra- 
dio sources (Jackson et al. [2007; Sowards-Emmerd et al. 2005 ) 
and their variability is hence characterized by blazar activity 
(Sect. 14.4] ). Interestingly, both show (1) a trend of reddening 
when they become brighter and (2) a trend of increasing intrinsic 
variability (for definition see Sesar et al. 120071 ) with increasing 
wavelength. Such a behaviour is opposite to typical radio-quiet 
quasars and also to Sharov 21 (Sect. 13. It . 

Results from the Palomar-QUEST Survey were recently pre- 
sented by Bauer et al. (2009). 3 113 objects were identified in 
7 200 square degrees with fluctuation amplitudes > 0.4 mag on 
time scales up to ~ 3.5 yr. There are only a few objects show- 



ing maximum amplitudes > 2 mag up to 3.7 mag; all of them 
are blazars. Nearly all of the 14800 spectroscopically identified 
quasars in the data base have jumps < 1 mag; the highest value 
is 1.8 mag. 

4.2. A microlensing event? 

As the flare of Sharov 21 appears to be a singular feature in the 
long light curve, it is tempting to speculate that it originates from 
a rare event. Here we first discuss microlensing. 

Chang & Refsdal (1979) first suggested that the flux of a 
(macrolensed) quasar can be affected by a star crossing close to 
the line of sight with time scales of the order of a few months 
to several years. On the observational side, quasar microlensing 
was first identified by Irwin et al. (1989). Since then, consider- 
able progress has been made in microlensing simulations, and 
observations of significant microlensing have been reported in a 
number of systems (e.g., Pelt et al 1998; Koopmans et al. 2000; 
Chae et al. BUOTl Wisotzki et al. 2003, 2004, Cha rtas et al. l2004l 
Eigenbrod et al. l2006l Paraficz|2006l Sluse et al. 120071 




Fig. 8. The crowded field towards Sharov 21 on the 10"x 
10"cutout from the 4-m KPNO image in the I band (N up, E 
left). 

Sharov 21 is seen through M31, the high star density close 
to the line of sight is illustrated by Fig. [8] The quasar is the 
brightest, slightly elongated object in the centre, all other ob- 
jects are most likely stars in M 3 1 . Note also that the quasar ap- 
pears slightly elongated which points towards an object at a dis- 
tance S 0'.'5. Unfortunately, there are no archival Hubble Space 
Telescope observations of the field. Extensive imaging was per- 
formed with the WFPC2 in 2008 to create an accurate map of 
M31 microlensing (PI: A. Crotts) where Sharov 21 is, however, 
several arcseconds out of the field. A deep image of the field was 
taken with Subprime-Cam at the 8 m Subaru telescope in 2004, 
but the quasar lies exactly in the gap between the two adjacent 
fields 6 and 7. 
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From the M 3 1 mass model (Geehan et al. 2006 ) we estimate 
a stellar column density of E» ~ 170 M Q pc~ 2 towards the quasar. 
This value is ~ 2.5 10~ 4 smaller than the critical surface density 
Serf, = Dc 2 l(4nG) with D = D S /(D L D LS ) and D LS = D s - D L , 
where Ds and Dl are the angular distances of the source and the 
lense, respectively, G is the gravitational constant. The optical 
depth, i.e. the probability for the quasar to fall into the Einstein 
radius of a star in M 31, is t = 2*/£ CI j t ~ 2.5 10~ 4 . 

In the case of quasar microlensing by stars in a foreground 
galaxy with high optical depth the lenses do not act individually 
and the light curve is complex. For low optical depth (r ^ 0.5), 
however, microlensing can be studied in the single-star approxi- 
mation (Paczynski 1986 ), which is applied here. More precisely, 
the assumptions are made that both the lens and the source are 
point-like and that the relative motion of the lens is linear. Then 
the light curve is given by the magnification 

F Y Mt) u(t) 2 + 2 



F v,gs u{t) yju{t) 2 + 4 



(1) 



where u(i) is the angular distance between source and lens in 
units of the Einstein angle ©e, and F VjO bs(0 and F v gs are the ob- 
served monochromatic flux density at time t and the mean flux 
density in the ground state, respectively, in the B band. As Ds 
and Dl are given, the light curve depends only on the mass 
Ml, the relative transverse velocity v t of the lens, and the im- 
pact parameter, i.e., the minimum distance M m ; n between lens 
and source. The light curve of a high magnification event can 
be significantly modified by the finite size of the source. This, 
however, occurs for u m [ n ~ R*/Re, with Re = ©eDl, whereas 
we have u m \ n ~ 10 3 R*/Re- Furthermore, for a scale of 8.4 kpc/ " 
at the redshift of Sharov21 and assuming that the source size of 
the UV radiation of the quasar is ^ 10 13 m, the source has an an- 
gular diameter about two orders of magnitude smaller than the 
minimum impact parameter and can be considered as point-like. 

The transverse velocity is determined by the motion of the 
lens in M31, the proper motion of M31 with respect to the 
barycentre of the Local Group (LG), the motion of the Sun 
around the Galactic centre, the motion of the Galaxy around the 
Local Group barycentre, and the motion of the LG relative to 
the cosmic microwave background (CMB). Since the first two 
effects are poorly constrained, we consider here for simplicity 
only the velocity of the LG with respect to the CMB. With 
vlg-cmb = 61 2 kms" 1 towards (l,b) = (270°, 29°) (Loeb & 
Narayan 2008) we have v t ~ 300 km s . 

For the simplifying assumption that all lenses have the same 
mass and velocity, the Einstein radius crossing time f£ = ©eM is 
constant and the event rate for quasar microlensing is estimated 
by F ~ 2 N q T/(ntE) (e.g., Mao l2008l ). where N q is the number of 
quasars. Given ~ 20 quasars with B < 20 per square degree and 
a surface area of ~ 2 square degrees within the 25 mag isophote 
of M 3 1 , we have N q ~ 40 and F ~ 1 per century. Hence, the 
discovery of a microlensed, faint quasar behind M 3 1 over an 
interval of half a century is not unlikely. 

When the tranverse velocity and the distances are fixed, the 
maximum amplification is determined by the impact parameter 
M m ; n and the time scale ?e is determined by M». For v t = 300 
km s _1 the light curve is best fitted with M„ ~ 0.3 M G and u m i n ~ 
0.048 (Fig.|9j top). A higher velocity requires a higher stellar 
mass. The fit is not perfect because the observed light curve 
shows a weak but clearly indicated asymmetry (Fig.|2| which 
hints on deviations from the single-lens hypothesis. Therefore 
we add a second lens which moves with the same v t but crosses 
the line of sight ~ 180 days earlier. For simplicity, each com- 
ponent is treated as a single star, the light curve is computed 
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Fig. 9. Light curve fi(t) of Sharov 2 1 around its maximum (no er- 
ror bars for lucidity, open symbols: uncertain data). Overplotted 
are the fitting curves (solid) for microlensing by one single star 
(top) and two single stars of the masses Me,\ and Ml,2 (bottom). 



via the multiplicative magnification approximation (Vietri & 
Ostriker 1983) pi — \i\ ■ fi2- A good fit to the shoulder in the 
light curve between 400 and 50 days before the maximum is 
achieved for (M %u M,, 2 ) = (0.3 M , 0.1 M Q ) if v, = 300 km s" 1 
and (1.2 M , 0.4 M ) if v t = 600 kms _I , respectively, with 
(«min,i, «min,2 = 0.055,0.8) in both cases (Fig.|9]i- The projected 
linear separation of the two lenses is ~ 35 AU. Both the mass 
ratio and the separation are not atypical of a binary star. 

Although microlensing by stars in M 3 1 appears to be a plau- 
sible explanation for the flare of Sharov 21, there are serious ob- 
jections. First, the probability for magnification as strong as in 
Sharov 21 is very low. With regard to his model with r = 0.1, 
Paczynski (1986) points out that "appreciable increases in inten- 
sity are there once per millennium ... there is not much hope in 
detecting intensity changes due to microlensing at such low op- 
tical depth". Second, colour index variations are expected in the 
case of quasar microlensing (Wambsganss & Paczynski 1991; 
Yonehara et al. 120081 ) if the source appears extended. For the 
point source-point lens constellation, however, an achromatic 
light curve is expected. Here, the source can be considered point- 
like, but B-R was considerably smaller in the flare compared to 
the faint state (Fig. [3). Third, as the angular separation of the two 
lenses is ~ ©e,i, the magnification pattern is expected to be more 
complex then for the single-star approximation. Detailed mod- 
elling with a more consistent treatment of the binary lens prob- 
lem is clearly necessary but is beyond the scope of the present 
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paper. It will be particularly interesting to see if such models 
find a natural interpretation for the steep rise of the flux before 
the peak. 

4.3. A stellar tidal disruption event? 

An alternative process which is rare and produces a strong 
UV/X-ray flare is the disruption of a star by the strong tidal 
forces of a massive black hole (Lidskii & Ozernoi 119791 Rees 
1988 11990b . Tidal disruption events (TDEs) were discussed so 
far mainly in the context of dormant black holes in non-AGN 
galaxies or in low-luminosity AGNs (e.g., Phinney 1989; Rees 
1990) and, more recently, of recoiling black holes (Komossa & 
Meritt 2008 ). At least a fraction of low-luminosity AGNs ap- 
pears to be powered by stellar tidal disruptions (Komossa et al. 
2004, Milosavljevic et al. 120061 ). For high-luminosity AGNs the 
situation is more complicated. TDEs are expected to be rare for 
black hole masses higher than a critical mass Mbh > M cl -i t « 
10 8 M Q where the gravitation (Schwarzschild) radius, R$, ex- 
ceeds the tidal disruption radius of solar mass stars so that such 
stars are swallowed whole without disruption (Hills 1975; see 
also Chen et al. [20081 Gezari et al. 120081 . On the other hand, 
it seems possible that a massive, self-gravitating accretion disk 
brings more stars into loss-cone orbits and enhances therefore 
the tidal disruption rate (Syer et al. |199U Donley et al. 120021 ). 

Until now, about a dozen TDE candidates in non-AGN 
galaxies have been found from X-ray surveys (Komossa 2002 ; 
Komossa et al. |2009l Esauei 120071 12U081 C appellu ti et al. 120091 
and also in the UV/Optical (Renzini et al. 119951 Gezari et al. 
2006, 2008 ). None of the events detected so far were found to 
be related to a high-luminosity AGN. Such a detection would 
be interesting because it provides an opportunity to check the 
basic tidal disruption theory as their predictions generally de- 
pend on Mbh which can be estimated independently in this case. 
However, as noted by Gezari et al. (120081) . the existence of var- 
ious mechanisms for the UV variability of quasars makes the 
interpretation of a UV flare subject to careful analysis. In partic- 
ular it is necessary to make sure that the observed flare is not just 
a more or less usual feature in a strongly variable light curve. 

The theoretical frame for the interpretation of TDEs has been 
set with the pioneering work by Hills ([1975), Lacy et al. 0982), 
Rees (fT988l . Phinney et al. ( [T9891 . Evans & Kochanek ([19891 
and more recently by, among others, Magorrian & Tremaine 
([19991 , Ulmer (fT^99b . Aya l et a l. (120001 Ivano v & Novikov 
d200"TTl, Men ou & Quataert (120011), L i et al. (120021, Bog danovic 
et al. ( 120041 , Wang & Merritt d2004l . Chen et al. ( 150081 . Lodato 
et al. ( 120081 ). We start with a short summary of the basic theory. 
A star of mass M* and radius originally in hydrostatic equi- 
librium and on a parabolic orbit, passing at the time fd within the 
tidal disruption radius 



Rl = R,(M hh /M t ) 



1/3 



(2) 



of the black hole will be torn apart by the tidal forces. The 
change in the black hole potential over the star will produce a 
spread in the specific energy of the gas. After the encounter, the 
gas particles will have a considerably widened energy distribu- 
tion. Following Rees d 1988b . about half of the mass has negative 
energy and is thus gravitationally bound to the black hole, the 
other half is unbound and will be ejected from the system. The 
bound elements move in highly excentric elliptical orbits of the 
size 2a = R p + /? a ~ R a where R p and R d are the pericentric and 
the apocentric distance, respectively. With the Keplerian relation 
between the energy and the size of the orbit and after replacing 



a by the orbital period P we have E = -\(2nGM bh /P) 2/3 for 
Mbh s> M». After one orbit the bound elements will return to 
R p and efficiently loose their energy and angular momentum via 
stream-stream collision and finally be accreted to the black hole. 
The consequence is a luminous flare at fn with a peak in the 
UV/X-rays and a spectrum which is characterized by a black- 
body temperature of a thick disk or a spherical envelope at the 
tidal radius R t from Eq. (0: 



T e <r — 



( L ) 


1/4 




-1/6 
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This scenario implicates two important conclusions for the 
light curve. First, if the energy distribution is uniform the mass 
return rate is determined by the relation between energy and pe- 
riod of the orbit dM/dP oc dE/dP oc P~ 5 / 3 . Assuming the time 
scale of the transformation of the orbital energy into radiation 
is short, i.e. the luminosity of the flare follows the accretion rate 
and the latter is given by the mass distribution of the return times, 
the flux density in the flare should be 



F oc At 



-5/3 



(4) 



where Af = t - t& is the time since the first passage of the peri- 
centre. Numerical simulations have shown that this 'standard' 
Af -5 / 3 light curve is a good approximation, at least for later 
stages (Evans & Kochanek 1989). Close to the peak luminos- 
ity the light curve can be substantially shallower (Lodato et al. 
120081 . 

As another consequence, the time Afo = fo - h the most- 
tightly bound material needs to fall back to R p is directly related 
to the mass of the black hole 



Afn 



10- 4 yr 



M bh R 

Mlk 3 



3\ 



1/2 



P = 



R 



p 



(5) 



with Mbh, M*,R* in solar units; k depends on the spin-up state of 
the star with k ~ 3 for the likely case that the star is spun up to 
near break-up spin and k ~ 1 if spin-up is negligible. 

In what follows, we shall check whether the total energy in 
the outburst of the Sharov21 light curve, the decline of the out- 
burst, and the time scales are consistent with the TDE scenario. 
The energy release in the flare is related to the mass of the star, 



-Eflare = ^/accM.C , 



(6) 



where r\ is the efficiency of converting mass to radiated energy 
and / acc is the fraction of mass of the star accreted to the black 
hole. For simplicity we assume r]f. dcc — 0.1, which might be 
quite high (see Li et al. 2002 ) but is not implausible, namely for 
/ acc ~ 0.5 following Rees ( 119881 ) and t] ~ 0.1 ...0.4 depending 
on the spin of the black hole. E^ m is obtained by integrating the 
bolometric luminosity over the flare where Lb o i(0 is computed 
from the monochromatic flux density, F v ^ ws ,(t) = F v,obs(t)-F v ^, 
with F v ,obs(0 an d ^v,gs in the restframe. A black body spectrum 
is assumed with T e g as a free parameter. The lowest possible 
stellar mass is M* ~ 6M G (.Eflare ~ 2 10 54 erg) corresponding to 
T e g ~ 2 10 4 K where the black body spectrum peaks in the B 
band, i.e., at A = 4400 A/(l + z) restframe. For such a spectrum 
a colour index B - R ~ 0.5 mag is expected after foreground 
reddening, which is in line with B - R observed for Sharov21 in 
the flare (Fig. [3). 

With r e ff , Lboi, and M, fixed, the stellar radius is constrained 
by Eq. (0 for a given black hole mass Mbh- In Fig. [10] (top), the 
resulting (M*,R,) combinations are compared with the stellar 
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Fig. 10. Stellar radius R t /R (top), ratio of tidal radius and grav- 
itational radius R t /Rs (middle), and return time Afo in days (bot- 
tom) as a function of the mass M* of the tidally disrupted star 
for M b h = 1, 2, 5, 10 1O 8 M (top to bottom). In the upper panel, 
the curves increasing from left to right are the R*-M t relations 
for the stellar models from Salasnich et al. (2000) for the base 
(bottom) and the top (top) of the red giant branch, respectively, 
for solar metallicity (dotted: solar abundance ratios, solid: a- 
enhanced abundances). The dashed horizontal line in the bottom 
panel marks Afo from the observed light curve (Fig.fTTI). 



models from Salasnich et al. (2000). For M bh ~ (1 . . . 10) 10 8 M 
we find M« ~ (10 ± 3)M G and R* ~ (200 ± 100) R @ which 
excludes main-sequence stars but not giants. For such giants, the 
tidal radius is clearly out of the gravitation radius, R$, of the 
black hole (Fig. [10] middle), whereas main-sequence stars of this 
mass range are again excluded as R t < R$. 

In order to check whether the decline from the maximum 
follows the 'standard' Af~ 5 ' 3 -law we perform linear regressions 
of log Fyflmre as a function of log At (Fig-fTTTl. Whereas the onset 
of the flare at fo is defined by the strong increase in the light 
curve at JD ~ 2 448 897, the beginning of the tidal disruption, 
ti, has to be considered as a free parameter. The resulting slope 
y = dlogF v fl are /dlogAf depends on A fo with y ~ -1.53-0.1Afo 
(days). The observed light curve is reasonably fitted by a power 
law with y — -5/3 for Afo, bs ~ 1.3 days and is shallower at the 
beginning of the burst, as predicted by Lodato et al. (2008). 

The bottom panel of Fig. [10] shows the fallback time Afo 
as a function of M, for different values of the penetration fac- 
tor B. Note that Afo from Eq. ([5]) is independent of Mbh be- 
cause Rl oz MZ} for given (T e ff, L, M t ) (Eq.[3]l. The short ob- 
served return time requires /3 ~ 30 for M, = 1OM . With 




1 1.2 1.4 1.6 

lop CAt/dai_|s5 




50 100 

At (days) 



Fig. 11. Top: Linear regression of the decline from the maxi- 
mum in the double-logarithmic presentation of the light curve in 
the quasar restframe for the Afo yielding a slope of y = -5/3. 
Bottom: Fit of the light curve in linear presentation by the result 
from the regression in the top panel (open symbols: uncertain 
data; no error bars for lucidity). The vertical dotted line marks 
the onset of the flare at fo. 



R p /R s = (4.8,2.4, 1.0, 0.5) for M bh = (1, 2, 5, 10) 1O 8 M and 
M, ~ 8...12M , the fundamental condition R p > R$ is 
matched for M b h < 5 1O 8 M . Guillochon et al. (2009) present 
highly resolved three-dimensional simulations of a tidally dis- 
rupted 1 M solar-type star approaching a 1O 6 M black hole 
with B — 7, i.e. R p /Rs ~ 3. Their results deviate from the case 
B — 1 used in previous simulations. In particular, the simulated 
light curve shows a phase of smooth increase before the abrupt 
rise to the maximum, which is qualitatively in agreement with 
the behaviour of Sharov 21 (Fig.[2]i. 

The TDE scenario for the flare of Sharov 21 is supported 
by (1) total energy in the outburst, (2) the shape of the light 
curve, and (3) the time scale of the event. However, it must be 
noticed that the standard TDE theory refers to inactive black 
holes. The presence of a massive, self-gravitating accretion disk 
is expected to have a significant effect on the dynamics of the 
stellar tidal debris. Moreover, the description of the stellar orbit 
for an encounter with the pericentre at only a few Rs requires 
a fully general-relativistic treatment which is, however, beyond 
the scope of this paper. 

One may argue that the probability for the capture of a 
~ 1OM giant star by a supermassive black hole is very low. 
Though this is certainly true, we emphasize that Sharov 21 -like 
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quasar flares are obviously extremely rare. Moreover it is worth 
mentioning (1) the high star formation rates and (2) the high 
fraction of massive stars in quasar host galaxies. It has been 
suggested for a long time that a luminous AGN phase is ac- 
companied with or follows after an intense starburst, especially 
when both kinds of activity are induced by a galaxy merger (e.g., 
Sanders et al. [19881 Granato et al. 120041 Springel et al. [2005 ). 
Lutz et al. d20081 > derived strong evidence of the presence of in- 
tense star formation (up to ~ 3 000 M & yr ) from Spitzer PAH 
detections in 12 type 1 quasar host galaxies at z ~ 2. Also from 
Spitzer IRS observations, Shi et al. (2009 ) derived star formation 
rates on the level of luminous infrared galaxies (10 . . . 1OOM 
yr ) for a complete sample of 57 type-1 quasar hosts at z ~ 1. 
This high star formation activity appears to be concentrated in 
the circumnuclear region. In nearby Seyfert galaxies nuclear 
starbursts have been observed within < lOOpc from the AGN 
(e.g., Davies et al. [20071 Watabe et al. 120081 . It is likely that 
more luminous AGNs are accompanied by more luminous star- 
bursts (Kawakatu & Wada 2008). The 'radius of influence' of the 
black hole is R\ = GM^/crl, where cr„ is the velocity dispersion 
of the surrounding stellar population. The Mbh-o"» relation from 
Merritt & Ferrarese d200TT > yields R m ~ 12(M bh /M ) 58 ~ 30 pc 
for M b h ~ 5 1O 8 M . The typical distance a 1OM giant can 
travel with a velocity v = cr» during its lifetime is ~ 1 kpc. (Note 
that the M^-cr, relationship is not a strong function of redshift 
up to z ~ 3: Shields et al. 120031 ) 

Moreover, the central region of a galaxy provides a peculiar 
environment for star formation as is indicated by the high num- 
ber of supergiants in our Galactic centre (e.g., Krabbe et al. l 199 1 1 
Najarro et al. 119941 Martins et al. 120071 Mauerhahn et al. 120071 1 
with nearly a hundred massive stars within the central parsec. 
Strong tidal forces, mass segregation, and other peculiar condi- 
tions may result in a relatively high Jeans mass for collapsing 
cloud cores (Morris & Serabyn 1996) and in a top-heavy initial 
mass function (Figer et al. 120021 Paumard et al. 2006 1 Bartko et 
al. 120091 . 

Gopal-Krishna et al. (12008 ) argue that the gas from TDEs 
may be a major factor for the abortion of radio jets in quasars. If 
the radio loudness is related to Mbh, as seems to be indicated by 
observations, the radio-loudness dichotomy can be explained by 
the existence of a critical black hole mass for the tidal disruption 
of solar type stars. As noted by these authors, the tidal disruption 
of giant stars is then required to explain the association of a few 
radio-quiet quasars with black holes of high masses. 

4.4. Other interpretations 

In the pioneering work of Cannon et al. ( 1 19681 119711 ), the class 
of quasars with variability amplitudes of the order of one mag- 
nitude and with often rapid fluctuations were called optically 
violently variables (OVVs). OVV AGNs are all radio-loud and 
their strong variability is believed to be due to relativistic beam- 
ing (Gaskell & Klimek l2003l and references therein). According 
to the current zoology of AGN types, OVVs constitute a sub- 
type of blazars, i.e. of jet-dominated active galaxies viewed close 
(~ 15°) to the axis of a relativistic jet (Urry & Padovani 1995). 
Viewing the jet almost directly head-on results in a great mag- 
nification of variations in the flux from the jet. Blazars vary dra- 
matically on a wide range of time scales from hours to years. 
In this context, the colour dependence of the variability found 
for the two strongest variable OVV quasars from the SDSS S82 
(Sect. l4.ll ) can be easily understood assuming that the spectral 
energy distribution is dominated by the synchrotron peak in the 
radio and infrared. This is the case for low-peaked BL Lac ob- 



jects and is typical also for flat-spectrum radio quasars: When 
the synchrotron emission from the jet makes a significant con- 
tribution which increases with wavelength, the optical spectrum 
becomes redder when the jet is brighter and the variability is 
stronger at longer wavelengths. Here we argue that Sharov21 is 
not an OVV quasar in the 'classical' sense: (1) The strong flux 
variation is limited to one single event, there is no evidence of 
strong variability on a wide range of time scales. (2) The spec- 
tral behaviour of the variability is typical of radio-quiet quasars 
(Sect.l4~TTi. (3) Sharov21 is not radio-loud. 

Variability of non-blazar AGNs has been attributed to 
various processes like oscillations in the accretion flow 
(Igumenshchev & Abramovicz 1999 ) or in the jet (Hughes et al. 
1998). Both models produce quasi-periodic features in the light 
curve where the maxima have approximately the same widths 
as the intervals in between, which is contrary to what we see in 
Sharov 2 1 . Avalanche flows in the accretion disk provide another 
explanation (Tekeuchi et al. 1995). The profile of a single flare in 
simulated light curves (Kawaguchi et al. [1998 ) is characterized 
by a gradual increase in brightness followed by a sudden decline. 
It is not clear whether the rapid increase to the maximum can be 
explained by such a model. 

An alternative explanation of AGN variability is based on the 
superposition of uncorrelated events such as supernova explo- 
sions (Terlevich ll992l) . star collisions (Courvoisier et al. 1996; 
Torricelli-Ciamponi et al. 2000), or colliding clumps in a clumpy 
accretion flow (Courvoisier & Tiirler 2005). First, we compare 
the energy in the flare of Sharov 21 with the radiated energy of 
supernovae. SN 2006gy has been considered the most luminous 
supernova with a total radiated energy £sN,rad ~ 10 51 erg (Smith 
et al. 120071 Ofek et al. 120071 ) . For the flare of Sharov 21 we esti- 
mate fiflare ~ 10 54 55 erg, dependent on the spectral energy dis- 
tribution, hence Efi. dre ~ lO^^ZssN.rad.max- It was the aim of the 
collision scenario to explain both the ultraviolet emission and 
its variability. As a modified version, such events could be rare 
and contribute flares to a continuum source the emission and 
variability of which are dominated by other processes. In fact, 
the star collision scenario is distantly related to the TDE model. 
The maximum amplitudes of the flares in the light curves pre- 
sented by Torricelli-Ciamponi et al. are much too low compared 
to Sharov21. However it remains to be tested if the model pa- 
rameters can be adjusted to fit our light curve. 

In another class of models strong flares are explained by pro- 
cesses in binary black holes. Lehto & Valtonen (1996) suggest 
a model where they associate the quasi-periodic outbursts of the 
OVV quasar OJ 287 to the times when the secondary crosses the 
accretion disk of the primary black hole. When this happens the 
secondary pierces a channel through the disk, the gas is heated 
in the channel, flows out, and radiates strongly. The model was 
demonstrated to be successful for OJ 287 (Valtonen et al. 2008 ). 
Other models attribute radiative outbursts to the tidal perturba- 
tion of the accretion flow due to the secondary (Sillanpaa et al. 
1988), or the sweeping of a processing relativistic beam across 
the fin e of sight, in analogy to SS 433 (Katz [19971 Liu & Chen 
120071) . OJ 287 is the best-studied quasar with a quasi-periodicity 
in the light curve. Its outbursts come in 12 yr intervals corre- 
sponding to 9 yr in the restframe. Our light curve covers, in the 
restframe, only ~ 10 yr before the outburst and ~ 5 yr after. If 
orbital time scales longer than about one decade are not atypical 
of supermassive black hole binaries, it cannot be excluded that 
the flare of Sharov 21 is related to such a process. It would be 
important to find more historical observations in plate archives 
which yield an extension of the light curve further back. This 
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however is much more difficult than for OJ 287 which is ~ 5 
mag brighter than Sharov 2 1 . 

Finally, we briefly mention the possibility that the quasar was 
by chance superposed by a nova in M 3 1 . This hypothesis was 
tested already by Sharov et al. ( 119981 ). From the measurement of 
the position before, during, and after the flare they find that the 
derived coordinates are in excellent agreement with each other 
and conclude that the flaring object is "coincident with the star 
that was visible at its place in 1982-1997". 

5. Summary and conclusions 

Sharov 2 1 , originally discovered by Nedialkov et al. ( 1 19961 ), was 
previously classified as a remarkable nova in M 3 1 (Sharov et 
al. I19981 ). Here we have shown that it is a typical type 1 quasar 
at z = 2.109 seen through the disk of M 31 and showing an ex- 
ceptionally strong UV flare. We created a significantly improved 
long-term light curve based on archival data, data from the liter- 
ature, and targetted new observations. Compared to the original 
data given by Sharov et al., the new light curve has a ~ 20 yr 
longer baseline and a better sampling. Altogether, more than 10 3 
single exposures from 15 wide field telescopes are included re- 
sulting in detections at 221 epochs from 1948 to 2009, with a 
relatively good time coverage after 1961. The data material is 
completed by a large number of observations without detection 
of the quasar but with useful upper brightness limits for the time 
interval from 1900 to 2009. Based on this data, we subdivide 
the light curve in two phases: a rather quiet ground state with 
B ~ 20.5 for at least 98% of the time and a strong outburst in 
1992 with an increase of the UV flux by a factor ~ 20. The vari- 
ability in the ground state does not significantly differ from that 
of other radio-quiet quasars of comparable redshift and lumi- 
nosity. A black hole mass of 5 10 8 M is estimated from the C iv 
line, corresponding to an Eddington ratio of ~ 0.6 for the ground 
state. By the comparison with ~ 8 000 quasars in the stripe 82 
of the SDSS on a 7 yr-baseline (Bramich et al. 2008 ) and more 
than 300 VPMS quasars (Meusinger et al. [20021 120031 with 
light curves having a time-baseline comparable to Sharov 21 we 
have demonstrated that the strong UV flare of Sharov 21 is very 
unusual for radio-quiet high-redshift quasars. We conclude that 
such a rare feature is the result of a rare event. As such we sug- 
gest two scenarios: (1) gravitational microlensing due to a star 
in M 3 1 and (2) a tidal disruption event (TDE) of a star close to 
the supermassive black hole of the quasar. 

In the TDE scenario, the total energy in the outburst can be 
explained by the disruption of a ~ 1OM giant star if we 'op- 
timistically' assume that half of the disrupted star is accreted 
to the black hole and ~ 20% of the accreted energy can be ra- 
diated away. The flare profile shows a sudden increase to the 
maximum followed by a decline which is reasonably fitted by 
the Ar 5 ^ 3 power law predicted by the standard TDE model. 
The short time span between the beginning of the tidal disrup- 
tion and the beginning of the flare, Afo ~ 1.3 days, requires an 
ultra-close encounter with B ~ 30 corresponding to a stellar or- 
bit with the pericentre at only a few R$- The present study has 
not taken general-relativistic effects into account. Moreover, it 
is unclear how the TDE scenario is affected by the presence of 
a massive accretion disk. Microlensing by a star in M31 is a 
plausible alternative explanation. Though the detection of a mi- 
crolensed quasar behind M 3 1 over half a century is not unlikely, 
high-amplification events corresponding to the flare of Sharov 2 1 
are very rare. We apply the point source-point lens approxima- 
tion to model the light curve. Assuming a transverse velocity of 
300kms _I , an acceptable fit is achieved for a low-mass binary 



with 0.3 M and 0.1 M Q . The observed light curve of Sharov 21 
is roughly fitted by either of the two scenarios, but more detailed 
modelling is necessary to decide if the flare can be reproduced 
accurately. Finally, we cannot exclude that the flare is part of a 
quasi-periodic activity similar to OJ 287 (Valtonen et al. 2008 ) 
on an intrinsic time scale of ^ 10 yr. The remarkable quasar 
Sharov 21 obviously merits further efforts, both for the comple- 
tion of the light curve and for its modelling. 

Finally, we notice that there is an interesting application 
of Sharov 21 -like flares. Projects like PTF, Pan-STARRS, and 
LSST will probably discover several such events. Due to their 
enormous luminosity and long time scale in the observer frame, 
such flares can provide background light sources for interven- 
ing matter and create thus interesting opportunities for high- 
resolution spectroscopy of matter at large distances. This argu- 
ment has been presented by Quimby et al. (2009) in the context 
of the brightest supernovae from the PTF but holds even more 
for quasar-flares comparable to that one of Sharov 21. 
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